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We obtain explicit susceptibilities for the countable infinity of phase transition temperatures of 
Muller-Hartmann-Zitartz on a Cayley tree. The susceptibilities are a product of the zeroth spin 
with the sum of an appropriate set of averages of spins on the outermost layer of the tree. A clear 
physical understanding for these strange phase transitions emerges naturally. In the thermodynamic 
limit, the susceptibilities tend to zero above the transition and to infinity below it. 


INTRODUCTION 

Forty years ago, Miiller-Hartmann and Zitartz [H-I^ 
showed that the Ising ferromagnet on a Cayley tree has 
a countable infinity of phase transitions of an unusual 
type. Specifically, in the low temperature ordered phase, 
when one traverses vertically across the zero-field region 
in the magnetic field - temperature phase diagram, they 
showed that the order of the phase transition can be 
anything between 1 and oo depending on the value of 
the temperature at that point. Therefore, they dubbed 
them phase transitions of continuous order. By identify¬ 
ing the values of temperature at which the order acquires 
integer values, they also extracted a countable infinity 
of phase-transitions in the low-temperature phase. The 
MHZ work constitutes one of the classic examples of a 
strange phase transition. Most of their work is rather 
technical and mathematical though, and a transparent 
physical understanding of the meaning of these strange 
phase transitions would be useful. 

Although criticized as unphysical, the Cayley tree (or 
the close cousin called Bethe lattice) is a popular geo¬ 
metric structure on which innumerable studies continue 
to be carried out in current times (3-[lflj| . One reason is of 
course the availability of exact solutions for some models. 
Another important, not often emphasized reason is that 
it simultaneously contains one-dimension like properties 
and infinite-dimension like properties, thus providing an 
interesting framework for theoretical explorations. For 
example, the partition function of the ferromagnet on 
a Cayley tree [III [l2| is identical to that of the one¬ 
dimensional chain (barring an inconsequential constant 
factor), and yet a phase transition exists in the Cayley 
showing infinite-dimensional character. 

MHZ extract the order of the phase transition at var¬ 
ious points in the low-temperature phase by studying 
the leading-order singularities of the full free-energy as 
a function of field as it is taken to the zero limit. Here 
we show that a direct transparent understanding of these 
phase transitions maybe obtained by the consideration 
of the pure zero-field model, building on a ‘memory- 
approach’ that we emphasized in recent work 0 . We 
were also partially propelled by the recent exact solu¬ 
tion of a one-dimensional long-range ferromagnet which 



FIG. 1. Schematic of the countable infinity of phase transi¬ 
tions of MHZ. The dots correspond to the countable infinity 
of phase transitions (as one traverses vertically by varying H 
through 0 at a given temperature in the ordered phase) of 
integer order. 

admits an unusual phase transition of a different type, 
namely mixed-order, which can simultaneously show 
a discountinuous jump in maganetization (first-order 
like), and a diverging correlation-length (second-order 
like) [110. 

THE COUNTABLE INFINITY OF PHASE 
TRANSITIONS OF MHZ 

Fig [T] carries a schematic of the phase diagram of MHZ. 
The gist of their approach is to make a careful detailed 
study of the free energy F{H,T). In a generic ferro¬ 
magnetic system with a phase transition, F{H, T) is an¬ 
alytic at all points except in the region shown in red: 
[H = 0,0 < T < Tbp]. In this limit the free energy is 
given by [i| 

F{H,T) = F(0,r) + fregiH^,T)+A{T)\H\^,H ^ 0, 

( 1 ) 

where the regular part is a function of because of 
symmetry, and the leading singular part A(T) shows a 
power law behavior. The unusual aspect of the Cayley 
tree lies in the fact that the critical exponent k varies 
continuously from I at T = 0 to oo Q at the usual phase 
transition called the Bethe-Peierls transition Tbp- This 
behaviour is in sharp contrast to most commonly encoun¬ 
tered phase transitions where k remains a constant (= I 
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FIG. 2. Figure shows a Cayley tree of depth n — 3 and 
coordination number 2 = 3. 


for a first order transition, 2 for a second order transi¬ 
tion and so on). By studying the points at which k takes 
integer values, they identify a countable infinity of phase 
transitions which fit into the Ehrenfest classification of 
phase transitions of integer power-order. 

Here we point out that in fact a simple (albeit unusual) 
set of susceptibilities identify the countable infinity of 
phase transitions, and indeed a clear physical picture of 
the meaning behind the phase transitions comes natu¬ 
rally out of them. The susceptibilities turn out to be a 
product of the zeroth spin with the sum of averages of 
appropriately grouped spins in the outermost layer of the 
tree. Furthermore, we can work entirely with the zero- 
field model, with no requirement of complicated proce¬ 
dures or mathematical methods related to the application 
of a tiny field followed by taking the zero-field limit. 


THE MODEL 

Following the notation of a recent piece of work involv¬ 
ing the author [l^, we consider the following Hamilto¬ 
nian: 

'H = -J^ a^aj +HY^aj, ( 2 ) 

{*j> i 

where the sum involves pairs of spins which are adja¬ 
cent on the tree (Fig. [5]) with coordination number 2 and 
depth n. Ui are Ising variables which can take values 
±1, and J is taken to be positive to make it a ferromag- 
net. The solution, when the external field is = 0, is 
trivial and we quickly recall it. We introduce new bond 
variables = uiaj. The 9ij can take values ±1, which 
make them effectively spin variables too. Specifying all 
the Oij, and the spin at the root of the tree ctoj com¬ 
pletely defines the system. The Hamiltonian then takes 
the following simple form 

H = (3) 

With the problem now rehashed into one with non¬ 
interacting spins under the influence of an external 


magnetic field, the partition function is readily written 
down [l^ : 


Z{J)cavley = 2{2 COSh{P J))^” = 2(2 COSh(^ J))"^”! , 

( 4 ) 

where Ni, is the number of bonds and N is the num¬ 
ber of spins, and j3 = 1/T is the inverse temperature as 
usual in equilibrium statistical mechanics. It follows di¬ 
rectly 13.ll6-19j| that the correlation function between 


any two spins Ui, aj is given by: 

{(JiCTj) = tanh(/3J)‘^‘^ = 


( 5 ) 


where we define a = tanh(,5J) for convenience, and dij 
is the distance of the (unique) shortest path between the 
points i,j. 

When the external field H is non-zero, there is no sim¬ 
ple closed-form expression, however MHZ [l[ write down 
an infinite-series expansion for the free energy F{H,T), 
and by a careful, elaborate study of the order in field 
at which the leading singularity occurs as one crosses 
zero-field at low-temperature, they obtain the following 
countable infinity of transition temperatures: 


tanh 

Z = 1, 2, • • • , oo, with I being identified as the order of the 
phase transition within the Ehrenfest scheme. Ti = 0 has 
a first order phase transition and Too = Tbp with order 
infinity is the so-called ‘Bethe-Peierls’ phase transition 
which is the temperature at which the system first orders 
as it is cooled down from high temperature. Fig [3] shows 
a Monte Carlo simulation carried out on a finite-sized 
system of depth n = 8 and coordination number 2 = 
3. We studied the dependence of magnetization as one 
crosses from negative to positive magnetic field at some 
of the transition temperatures of MHZ. It seems plausible 
that in the thermodynamic limit, at a higher transition 
temperature, a greater derivative of the magnetization 
with respect to field would diverge at H = 0. Although 
the simulations were run on a finite-system the data at T 2 
display a considerably sharp drop near H = 0 indicative 
of the diverging derivative, since it is a second order phase 
transition. 

Here we show that the above countable infinity of 
transition temperatures may be directly obtained from 
the zero-field model bypassing the elaborate complicated 
methods of studying the infinite series and the order of 
divergence of the free energy in the limit oi H 0. We 
do this by an explicit construction of a special set of ‘sus¬ 
ceptibilities’ for the transitions Tj. 


CONSTRUCTION OF THE SUSCEPTIBILITIES 
FOR THE MHZ TRANSITION TEMPERATURES 

In order to construct the susceptibilities for the MHZ 
transition temperatures, we choose the depth of the lat- 
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FIG. 3. A plot of the magnetization m = versus 

magnetic field from Monte Carlo simuations of a finite-sized 
sytem of coordination number z = 3 and depth n = 8. Data 
are shown for different runs at the MHZ transition tempera¬ 
tures r 2 , 7io,Tiooo,Tho =Tbp. At higher transition temper¬ 
atures, the curve becomes smoother and smoother indicating 
that in the thermodynamic limit the singularity would occur 
at a higher order differentiation with respect to H. Since our 
system is finite, the data are practically indistinguishable for 
Tiooo and Too- 


tice to be of the form n = Ip, where I can take values 
1,2,3, ••• , oo. Since we are interested in the p —)• oo 
limit, no loss of generality is incurred. The number 
of spins in the layer of the tree is Af„ = z{z — 
1)71-1 = (1 _|_ i)y", where we have defined 7 = (z — 1 ) 

for convenience. Let us denote the layer spins by 
fn,!, 0 - 71 , 2 , • • • , 0 - 77 ,Ar„ in Order from left to right, as can be 
visualized in Fig. [2] Next, we group the first Mp = 7 ^* 
spins of the layer and call their average an,i = 
^n,i] we then group the second Mp spins of the 
n*'* layer and call their average an ,2 = Yh=m^+i 
and so on. By this procedure we form Kn,p = ^ spin 
averages: 

1 jMp 

= Jf H 

^ i=U-^)Mp + l 

J = 1 , 2 , 3, • • • , Kn,p- 

Now we are ready to write down the susceptibilities. 
Recalling that (Tq is the spin at the root node, the sus¬ 
ceptibilities are simply given by: 

Kp,p 

O"o ^ ^ ^n,j: (^) 

1 = 1 

which is our main result. To see that this leads to the 
MHZ transition temperatures, let us invoke the two-point 
correlation functions from the last section to compute the 


expectation value: 


Kr^ 


{Xi) = (cto X! 


i=i 


Kn,p - j^p 

-Si S 


((To 


Krt 


j^p 


"■ty -I ^ —y 

-Si s 


1=1 


i=0-l)Mp-|-l 


= i^77.pa" 

= (l + -)[7^'-^i'f. 
7 


(9) 


( 10 ) 


Therefore we see that as p —>■ 00 , 


(Xi) 



if yO-Hai < 1 
if yO-Hai > 1 . 


( 11 ) 


I , 

7 a= 1 thus defines the I transition temperature. 
These are precisely the transition temperatures of MHZ 
as given in Eqn. | 6 ] [l|, [ 2 ^ . It is worth pointing out that 
for Z = 2 , this yields the so-called (because of the ap¬ 
pearance of the same in the disordered version of the 
same Hamiltonian [l^ l ‘spin-glass’ transition tempera¬ 
ture Tsg = I— , _i 1 ' , and the I —>■ oo limit yields the 

Bethe-Peierls transition temperature Tbp = tanh~^(l-) ' 


CONCLUSIONS 

We have introduced a set of ‘susceptibilities’ that help 
identify the mysterious MHZ transition temperatures on 
the Cayley tree in a transparent manner. They are given 
by the product of the zeroth spin with an appropriately 
averaged sum of spins from the outermost layer in a Cay¬ 
ley tree. A clear physical understanding of the phase- 
transitions emerges naturally. We observe that our sus¬ 
ceptibilities have the feature that in the thermodynamic 
limit, they are zero above the phase transition, but tend 
to infinity below it. We are able also to identify the 
second-order phase transition T 2 as the phase transition 
known as the spin-glass phase transition in the literature. 

Furthermore, although we have concentrated on the 
Ising ferromagnet here, the susceptibilites defined here 
are primarily attached to the geometry of the lattice. 
Therefore, they should be applicable much more gener¬ 
ally: for example with m-component vector spins and 
the bonds could be ferromagnetic or antiferromagnetic 
or disordered. Quantum models should display similar 
transitions as well: a detailed investigation of various 
models from this perspective would be desirable. 
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Finally, we remark that the statistical mechanics prob¬ 
lem on the tree has been connected with the problem of 
reconstruction of information on trees in formal, exten¬ 
sive studies 21, . It would be interesting to under¬ 


stand if and how the MHZ countable infinity of phase 
transitions would fit into this generalized problem, which 
might be of interest to a broader community. 
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